Signatures of the Anderson localization excited by an optical frequency comb 
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We investigate Anderson localization of light as occurring in ultra-short excitations. A theory 
based on time dependent coupled-mode equations predicts universal features in the spectrum of the 
transmitted pulse. In particular, the process of strong localization of light is shown to correspond 
to the formation of peaks in both the amplitude and in the group delay of the transmitted pulse. 
Parallel ab-initio simulations made with Finite-Difference Time Domain codes and molecular dy- 
namics confirm theoretical predictions, while showing that there exist an optimal degree of disorder 
for the strong localization. 



Introduction — In the last years, progress in laser tech- 
nology has resulted in the generation of pulses with du- 
ration of few optical cycles [lH4j, whose spectral con- 
tent is an optical frequency comb, i.e., a discrete, regu- 
larly spaced, series of sharp lines in the region of opti- 
cal frequencies. These new laser sources have favoured 
many outcomes in fundamentals research and applica- 
tions, ranging from optical frequency metrology [5|-lZ| to 
pulse synthesis [1, @]. In the same years, the system- 
atic study of light propagation in disordered dielectrics 
has stirred particular attention In these systems 

the multiple scattering can become so strong that the 
diffusive approximation for the photon transport lapses 
and the randomness of the three-dimensional (3D) struc- 
ture leads to the emergence of Anderson localizations, 
i.e., long living exponentially localized electromagnetic 
resonances. In several respects, the observation of these 
localizations is still largely debated [18l-l24j. One of the 



leading issue is determining universal features of light lo- 
calization dynamics in order to infer their excitation from 
the analysis of the spectrum (amplitude and phase) of the 
transmitted pulse. In one-dimension (ID) transmission 
fluctuations can be related to exponentially localized or 
to delocalized necklace states, both mediating a resonant 



like behaviour (see, e.g., [25|, [26|). In 3D, with refer- 



ence to dielectric samples, emphasis has been given to 
the shape of the tail of the transmitted pulse, that should 
display deviations from a diffusive exponential trend in 
the presence of long living states [19j. However for a 
broad-band excitation such an analysis cannot be directly 
applicable because: (i) the diffusion approximation may 
not hold true in whole excitation spectrum, (ii) material 
dispersion should be accounted for. In addition, when 
dealing with ultra- fast laser sources, one could argue on 
which is the Frequency Resolved Optical Signal (FROG) 
in the presence of Anderson localization (see, e.g., |27j), 
or which is the corresponding amplitude/phase profile 
in the spectrum and if they display universal distinctive 
features. In this respect, a time-resolved electromagnetic 
approach is unavoidable, especially for dealing with a fre- 
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FIG. 1: (Color on line) Localization couplings and energy 
propagation scheme in the TDCMT. 



quency dependent refractive index and broad-band exci- 
tations. 

In this Manuscript, we investigate the response of a ran- 
dom medium to a frequency comb. We tackle this issue 
by Time Domain Coupled Mode Theory (TDCMT) [H|, 
and by an ab-initio parallel numerical approach that em- 
ploys finite- difference time-domain (FDTD) codes com- 
bined with a Molecular Dynamics (MD) [29|, which allow 
to completely simulate the dynamics of Maxwell's equa- 
tions (FDTD) inside realistic configurations of scatterers 
(MD). We compare materials with different index con- 
trast, and, for a strongly scattering medium, we vary 
the filling fraction as a measure of the degree of ran- 
domness. Theory and simulations demonstrate that the 
transmitted pulse contains spectral universal features of 
the Anderson localization of light, which, according to 
the FDTD simulations, is mostly evident at a specific 
amount of disorder. This is the first ever reported first- 
principle study of the role of the 3D Anderson localiza- 
tion of light in the presence of a flat wide-band excitation, 
and it is expected to stimulate further experimental and 
theoretical work. Indeed, so far, most of the reported in- 
vestigations dealt with a diffusive almost-monochromatic 
regime; and a full Maxwell solution for a frequency-comb 
input was not reported before. 

Theory — In our theoretical picture the impinging 
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laser pulse (propagation direction z, input surface at 
z = 0) excites a series of Anderson localizations at dif- 
ferent frequencies, which are expected in proximity of 
the surface of the sample; energy is then transported 
in the sample by coupling between localized states at 
the same frequencies. We focus on the output signal at 
the input propagation direction, thus only those localized 
modes that are overlapped with the beam contribute to 
the output, while the others behave as radiative losses. 
The coupling between resonances at different frequencies 
will not contribute to the energy transport, indeed their 
beating will average out. Such an approach can in prin- 
ciple applied to any dimensionality to describe transport 
of energy mediated by coupled localized resonances. We 
consider M chains of uncoupled oscillators, each corre- 
sponding to a resonance uj u (y — 1,2,...,M) with am- 
plitude a^ n with n = 1...N, where N is the number of 
Anderson localizations at frequency uj u overlapped with 
the input beam (FigJTJ [we group resonances at the same 
position z (within a localization length) into a single os- 
cillator a Uy7l \. The propagation of energy at frequency uo v 
is then given by the linear TDCMT equations: 



2 7^— = 

dt 

(1) 

with r n ^ accounting both for the emission losses in di- 
rections different from the input beam and for material 
absorption. The coeffcients fi in (pQ) describe the cou- 
pling between modes and, in general, they also account 
for material dispersion. For n = 1, the previous equa- 
tions contain /io,^o,iv = ^2/r e si(uj u )^ which is the frac- 
tion of the signal coupled at frequency uj u , being si(uj u ) 
the amplitude of the supercontinuum at uj u and 2/r e the 
coupling loss rate in z-direction (Fig. [T]) .28]. Since the 
input beam will exponentially decrease with the sample 
length, we assume that only one single mode for each uj v 
(the nearest to the input surface z = 0) will be directly 
excited by the input beam. For n > 1, a n -\^ v (a n+ i^) 
is the Anderson localization at uo v preceding (following) 
&v,n' The transmitted signal is then s-t = fijsj ydjs v , be- 
ing \i the coupling coefficient with the output region. The 
transmission T = St / Si is then found by solving (pp) in the 
Fourier domain (^ — » iuj): 



N 

k=i 



(2) 



being G(k,v) the following continued fraction: 



A*fc-i,i 



G(k,v) 



i(u-^)+^-+li k+ltt ,G(k+l) ' 

.0, 



k < N 
k > N 



(3) 



Universal properties of the system can be demonstrated 
from ([2])-(j3|), by looking at the two physical observables: 
the transmittance |T| 2 and the group delay ^ (with (j) 



being the phase of T). We begin by demonstrating with 
mathematical induction the property 9[G(ra, v)\ \ UJ = UJu = 
0. This is trivially fulfilled for G(N), and also for G(k): 

~^k-\,v [(u ~ u v ) + /ifc+i,„3?[G(fe + 1)]] 



9f[G(fc)] 



if 3f[G(Jfe 



•^) + ^ + M/c+i^G(A:+l)| 2 



(4) 



1,^ 



uj u = 0. In the same way we can 
verify that ^[ dG ^J u " > ]\ u=u;u = and, as a corollary, 



Z[G(k,v) 2 ]\ u 



0, U[G{k,vf d ^l]l 



and 



^ d^ 



UJu = 0. By induction we then obtain: 
d\G{k,v)\ 2 . „ d 2 ^ 



0, 



0, 



du> '™ ~' du 2 lu=UJ » (5) 

being <pk,v the phase of the continued fraction G{k,v). 
Correspondingly \T\ 2 ex Y\ k \G(k,v)\ 2 and (j) oc ^2 k (/>ky, 
hence, owing to Eq.(j5j), we are able to say that in the 
presence of Anderson localization occurring with multi- 
ple resonances, as produced, e.g., by broadband excita- 
tions, for each resonance uj v the transmission spectrum 
exhibit extrema in the amplitude and the group delay. 
Sample and numerical setup — The previous general ar- 
guments have a universal character and applies to a vari- 
ety of different situations; in order to make a specific 
example and verify the validity of our theoretical ap- 
proach we resorted to first-principle FDTD parallel sim- 
ulations. To consider a realistic sample of disordered col- 
loidal system we used MD simulations, which furnish a 
distribution of 8000 mono-dispersed spheres with radius 
r and filling fraction <£> in air interacting with a Lennard- 
Jones potential. We first consider two different materials 
with r ~ 150 nm at a filling fraction $ ~ 0.55: Silica 
Si02 glass (n = 1.5) and Titanium Dioxide Ti02 in its 
Rutile form (n = 2.5). With reference to realistic fre- 
quency comb generators, we tailored the temporal pro- 
file of the impinging pulse so that its spectral content 
is centred around the carrier frequency fo = 375 THz 
(Ao = 800 nm), and regularly spaced between 320 THz 
and 430 THz by a repetition rate f rep ~ 90 GHz. The 
input pulse E(t) = sin(27r/ n t), where the sum is per- 
formed over 1000 equi-spaced frequencies f n (Fig|2j). 

In our numerical experiments we take a train of ultra- 
short, y-polarized pulse impinging on the x-y face of the 
sample at normal incidence. The input spatial profile is 
Gaussian TEMoo with waist wo = 1 jiva. 10 ps runs take 
one hour by using 1024 processors on an IBM system. 

For each simulations, we collect: i) the spatial profiles 
of the electric field at the output plane of the sample; ii) 
the electromagnetic energy 6 and iii) the corresponding 
spectrum. Afterwards we calculate: a) the group delay as 
the derivative, with respect to frequency, of the phase of 
Fourier transform of the electric field; b) the frequency- 
resolved optical gate (FROG) signal Ifrog{u,t) as: 



lfrog 



(u,t) oc | 



E y {t)Ey{t-T)e- iu3t dt\ 2 (6) 
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FIG. 2: (Color online). Ti02 refractive index n vs frequency, 
as obtained from the Sellmeier equation after Devore [30! ]. 
The filled area corresponds to the frequency comb spectrum. 
The inset shows the input pulse. 
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FIG. 3: (Color on line) Slices (x, z) (in the middle of y axis) of 
the electromagnetic energy £, resulting from dispersive simu- 
lations of Si02 (a) and T1O2 (b). Corresponding FROG signal 
of the transmitted pulse (E y ) for SiCb (c) and for TiOi (d). 



where E y (t) is the output signal, which corresponds to 
the complex signal s t from TDCMT above. 

Numerical results — The first signature of the exis- 
tence of localized modes in Ti02 is found by looking at 
Fig|3^,b , where we notice a drastic change in the spatial 
profile of the Energy £ from extended states in Si02 ([3fci) 
to localized states in TiC^fO}). To characterize the pres- 
ence of localized modes, we compute the FROG signal 
Ifrog(oj,t) [23, [HI, [32|. Figure [3] shows I FRO g(oj,t) 



calculated by the Eq.(|6]): panel (c) corresponds to the 
field E y (t) coming out by Si02, while panel (d) refers to 
Ti02- The FROG signal is able to detect if the strong 
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FIG. 4: (Color on line) Squared modulus of the Fourier trans- 
form of E y (dashed line, left axis) and group delay (continuous 
line, right axis) for Si02 (a) and T1O2 (b). The vertical axis 
of the group delay is shifted such that the zero corresponds to 
the exit time of the low frequency side of the input spectrum. 



localization is set in the medium. Indeed, the pattern in 
the right panel (d) displays a spectrogram in which all 
the frequencies (reported in the y-axis) survive for longer 
time (~ 1.2 ps) with respect to the case of Silica glass 
in panel (c) (~ 0.2 ps) [Note that the FROG spectrum 
is doubled with respect to the field spectrum]. This ob- 
servation is in agreement with the emergence for Ti02 of 
several localized modes with a long lifetime, correspond- 
ing to amplitude peaks and longer delays. Comparison 
with our theoretical predictions can be tackled through 
the spectral analysis on the Fourier transform T(uo) of 
the output pulse E y (t): 



T(uj) = A(uj)exp[i^^L} = A^e^l 



(7) 



The derivative dip/duj is equal to the product between the 
inverse of group velocity v g = du/dk and the distance L 
covered by the incident pulse at the end of the sample, 
that is the group delay r g = L/v g . 

Figure 3] shows the squared modulus A 2 (uj) = \T(uj)\ 2 
and the group delay r g (uj): we compare the results for 
Si0 2 (Fig. gt) with those Ti0 2 (Figgjb). The transmis- 
sion gives for the transport mean free path I = 0.6 fim 
for Ti0 2 (kl = 5) and I = 10 fim for Si0 2 (kl = 80). 
Localization is accompanied in Ti02 by the formation of 
several peaks for A 2 (uj) and r g (uj) (absent for Si02) in 
agreement with theoretical predictions. Note that due to 
the multiple couplings (i.e., from a resonant state to an- 
other or mediated by back-coupling to the input beam), 
both minima and maxima for the delay and the ampli- 
tude are obtained (corresponding to mixed side-cou plin g 
and forward-coupling regimes for optical cavities [28|). 
Such a response is a universal feature that we found for 
any considered disorder realization. We also compared 
dispersive and non-dispersive materials, which resulted 
into quantitatively different distribution of peaks (not 
shown), while retaining for both cases the predicted fea- 
tures. 

The link between the amount of disorder and the on- 
set of Anderson localization is a key issue. Previous nu- 
merical investigations for 3D photonic crystals [23] have 



shown that the localization length is minimal for a spe- 
cific degree of disorder. Notwithstanding the fact that 
we are considering a completely disordered regime, well 
beyond the perturbative regime considered in [23], here 
one expects that such an optimal condition corresponds 
to the largest fluctuations of the group delay. Hence 
we compared various Ti02 structures exhibiting different 
filling fractions, as obtained by varying the radius of the 
(eventually overlapping, as in porous media, see insets in 
figure [5] below) spheres. The filling fraction $ is adopted 
here as a measure of the degree of disorder, indeed, for 
low (high) <I> the structure is mainly by air (TiC>2), and 
can be taken as an ordered, homogeneous-like, system: 
an optimal value for <I> must exist for Anderson mediated 
electromagnetic resonances. 

In figure [5] we show the spectra attained for three differ- 
ence values of the group delay are clearly more pro- 
nounced in the panel (b). In the bottom panel we show 
the maximum group delay versus the filling fraction, that 
display the greatest values in the range G [0.3, 0.5], cor- 
responding to optimal disorder for the localization. Note 
that a modulated amplitude is also present for large $ 
(right panel), but the attained delays are much smaller 
than for <I> = 0.4 (or vanishing as in figure H^i above); 
this is related to boundary reflections (for large <3> the 
whole structure is a cube resonator) and "sub-critical" 
(i.e. weakly) localized stated. Indeed, Anderson local- 
ization is signaled by the presence of large modulations 
in both the transmission and in the delay; at variance with 
the ID case (see, e.g., [25]), in 3D a modulated trasmis- 
sion does not necessarily correspond to an Anderson lo- 
calization. On the other hand, the narrow and sharp flat 
region for the delay versus $ in Fig|5]signals the existence 
of two phases (i.e., diffusive and localized); large delays 
for a specific disorder are directly related to the onset of 
the strong-localization, which also results in the FROG 
signal in FigJH We stress that the attained maximum 
delays for the various $ in the "optimal region" corre- 
spond to different resonances, thus denoting a largely 
frequency dependent regime, if compared to a shallow 
diffusive propagation. 

Conclusions — We have reported on the universal fea- 
tures of the Anderson localization of light in the presence 
of an ultra-wide band excitation. We have shown that 
strong localization phenomena can be discriminated by 
the appearance of a resonant like response in the ampli- 
tude and by the group delay versus frequency, which is 
maximum at a specific disorder that can be interpreted as 
the optimal configuration for trapping light in a random 
material. 
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FIG. 5: (Color on line). (Top panels) Spectral response 
(squared modulus, left scale, dashed line, and group delay, 
right scale, continuous line) for different filling fractions for 
Ti02- (Bottom) Maximum attained group delay versus fill- 
ing fraction; the insets shows the material distribution, cor- 
responding to the top panels. 
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